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Abstract 

We investigate the dual superconductivity hypothesis in pure SU(2) lattice gauge 
theory. We focus on the dual Meissner effect by analyzing the distribution of the 
color fields due to a static quark- antiquark pair. We find evidence of the dual 
Meissner effect both in the maximally Abelian gauge and without gauge fixing. 
We measure the London penetration length. Our results suggest that the London 
penetration length is a physical gauge-invariant quantity. We put out a simple 
relation between the penetration length and the square root of the string tension. We 
find that our extimation is quite close to the extrapolated continuum limit available 
in the literature. A remarkable consequence of our study is that an effective Abelian 
theory can account for the long range properties of the SU(2) confining vacuum. 
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I. INTRODUCTION 

Understanding the mechanism of quark confinement is a central problem in the high 
energy physics. This requires, among other things, to identify the dynamical variables 
which are relevant to the confinement. 

A satisfying solution would be to set up an approximate vacuum state which confines 
color charges. This way one could derive an effective action which describe the long- 
distance properties of QCD []T|. Even this incomplete program, however, mandates a 
non-perturbative approach. Fortunately we have at our disposal a framework in which we 
can do non-perturbative calculations, namely the lattice discretization of gauge theories. 
Since a typical Monte Carlo simulation generates vacuum configurations, one expects to 
gain information on the non-perturbative vacuum structure. 

However a guideless search into the numerical configurations generated during Monte 
Carlo runs is hopeless. In other words, we need some theoretical input which selects the 
dynamical variables relevant to the confinement. The situation looks similar to the theory 
of superconductivity. Indeed, it was the Cooper's observation that the Fermi surface 
is unstable with regard to the formation of bounded electron pairs which led Bardeen, 
Cooper, and Schrieffer to formulate the successful BCS superconductivity theory 0. 

An interesting possibility has been conjectured long time ago by G. 't Hooft and 
S. Mandelstam 0]. These authors proposed that the confining vacuum behaves as a 
coherent state of color magnetic monopoles. This is equivalent to say that the vacuum is a 
magnetic (dual) superconductor. This fascinating proposal offers a picture of confinement 
whose physics can be clearly extracted. As a matter of fact, the dual Meissner effect causes 
the formation of chromoelectric flux tubes between chromoelectric charges leading to a 
linear rising potential. It is worthwhile to discuss briefly the 't Hooft's proposal [^. 

Let us consider the non-Abelian gauge theory spontaneously broken via the Higgs 
mechanism. The Higgs flelds are in the adjoint representation. For concreteness we focus 
on the Georgi-Glashow model It is well known that the Georgi-Glashow model allows 
field configurations which correspond to magnetic monopoles [j^. Moreover, one readily 
finds that the monopole mass is given by 

M^o, = C^, (1.1) 
a 

where Myy is the mass of the charged vector boson, C a constant and a the fine struc- 
ture constant. The dual superconductor scenario is realized if these magnetic monopoles 
condense by means of the magnetic Higgs mechanism. This means that the monopoles 
become tachionic: 

< . (1.2) 

mon — \ J 

From Equation( |l.l|) we see that — > (if we kept a fixed). The fact that must 
go through zero suggests that the original Higgs field could be removed. Thus we are led 
to consider the pure gauge theory without elementary Higgs fields. The role of the scalar 
Higgs field is played by any operator which transforms in the adjoint representation of the 
gauge group. More precisely, after choosing an operator X{x) which transforms according 
to the adjoint representation, one fixes the gauge by diagonalizing X{x) at each point. 
This choice does not fix the gauge completely; it leaves as residual invariance group the 
maximally Abelian (Cartan) subgroup of the gauge group. Such a procedure is known 
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as Abelian projection 0. For instance, if the gauge group is SU(N), then after gauge 
fixing the residual invariance group is U(l)^~^. The world line of the monopoles can be 
identified as the lines where two eigenvalues of the operator X{x) are equal. The dual 
superconductor idea is realized if these Abelian monopoles condense. 

It is evident that the monopoles are dynamical; they will take part in the djTiamics of 
the system. As a consequence the problem of monopole condensation cannot be dealt with 
the perturbation theory. On the other hand, the Abelian projection can be implemented 
on the lattice . Thus one can analyze the dynamics of the Abelian projected gauge 
fields by means of Monte Carlo simulations. In the following we shall consider the pure 
SU(2) gauge theory. 

To perform the Abelian projection we make a choice for X{x). The simpler possibility 
is to consider a local quantity. For instance, we can use a plaquette with a definite 
orientation (field-strength gauge) or the Polyakov loop (Polyakov gauge). In these unitary 
gauges we implement the gauge fixing by means of the matrices V{x) which diagonalize 
X{x) at each lattice site: 

V{x)X{x)V\x) = diag[e"'^''\e-"'^''^ . (1.3) 

It is straightforward to check that the residual gauge invariance group is the U(l) group 
with transformations exp(i(T3^(x)). 

The Abelian projection of the gauge transformed links 

U^{x) = V{x)U^{x)V^{x + i^) (1.4) 

amounts to write 

U,ix) = W,{x)U^{x) (1.5) 

with 

(1.6) 



diag 



arg 
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:i.7) 



U^{x) is the Abehan projection of U^{x). 

A different class of gauge fixing has been proposed in the literature, namely the Abelian 
covariant gauge or maximally Abelian gauge 0. In the continuum the maximally Abelian 



gauge corresponds to impose the constraints |10 



x) = 



:i.8) 



A], ± lA^,, and Dn is the -covariant derivative. On the lattice the con- 



where 

straints ( |1.8| ) can be implemented like the Landau gauge 
corresponds on the lattice to diagonalize 



la, 12 



Indeed Equation 



X{x) = Y,{U,{x)a,Ul{x) 



+ UUx-fL)a;U^{x 



A)} 



:i.9) 
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To do this it is enough to maximize iteratively the quantity 

R = J2[^A(^)(rsUl{x)] , (1.10) 

X,fl 

where the f/^(a;)'s are the gauge transformed hnks ( [L.4| ). We thereby obtain the matrices 
V{x) and perform the Abehan projection of the hnks by Eqs. ( |1.4|) -( pr7D . 

From the above discussion it is evident that the monopole dynamic does depend on 
the choice of the operator needed to fix the gauge. On the other hand the confinement of 
color charges via monopole condensation can not depend on the gauge fixing. However, 
it is conceivable that the dual superconductor scenario could manifest with a judicious 
choice of X{x). This outcome could arise from a gauge fixing which freezes the degrees of 
freedom which are irrelevant to the confinement. We feel that the situation is similar to the 
time-honored BCS theory of superconductivity. Indeed in the BCS theory one deals with 
a reduced Hamiltonian which breaks the electromagnetic gauge invariance. Nevertheless, 
the reduced BCS Hamiltonian offered the correct explanation of the Meissner effect. As a 



matter of fact, it was showed [0 that the collective states which are essential to restore 
the gauge invariance do not contribute to the BCS calculation of the Meissner effect. In 
other words, the reduced BCS Hamiltonian, by retaining the degrees of freedom relevant 
to the superconductivity, gives a sensible answer even though it breaks the electromagnetic 
gauge invariance. 

Interesting enough, it turns out that, if one fixes the maximally Abelian gauge, the 
Abelian projected links seem to retain the informations relevant to the confinement 



Thus, it is important to deepen the study of the dynamics of the Abelian projected fields 
in that particular gauge fixing. 

The aim of the present paper is to analyze the finger-print of the dual superconduc- 
tor hypothesis, namely the Meissner effect. To this end, we analyze the distribution of 
the color field due to static quark-antiquark pair in SU(2) lattice gauge theory in the 
maximally Abelian gauge. Moreover we will study the gauge dependence of the London 
penetration length. A partial account of this paper has been published in Ref. [^|. 



The plan of the paper is as follows. In Sec. II we explore the field configurations 
produced by the quark-antiquark static pairs both in the case of Abelian projected links 
after maximally Abelian gauge has been fixed, and in the case of full SU(2) links. In 
Sec. HI we analyze the transverse distribution of the longitudinal chromoelectric field. In 
Sec. IV we investigate the relation between the penetration length and the string tension. 
Our conclusions are relegated in Sec. V. The Appendix comprises several technical details 
on the maximally Abelian gauge fixing. 



II. COLOR FIELDS 

In this Section we analyze the distribution of the color fields due to static quark-antiquark 
pairs. Following the authors of Ref. [|16|, we can measure the color fields by means of the 
correlation of a plaquette Up with a Wilson loop W. The plaquette is connected to the 
Wilson loop by a Schwinger line L (see Fig. [|). Moving the plaquette Up with respect to 
the Wilson loop one can scan the structure of the color fields. In a previous study [17| we 
found evidence of the dual Meissner effect in the maximally Abelian gauge. In particular 
we measured the penetration depth of the flux tube chromoelectric fleld. However in 
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Ref. [0 we employed rather small lattices [L = 12). In this work we extend our previous 
study in two directions. Firstly, we perform numerical simulations on lattice whose size 
ranges from L = 16 up to L = 24. In addition we investigate the gauge invariance of 
the penetration length. To do this we perform the numerical simulations both in the 
maximally Abelian gauge and without gauge fixing. 



A. SU(2) 

According to Ref. ||16|, one can explore the field configurations produced by the quark- 
antiquark pair by measuring the connected correlation function (Fig. |l|) 

{ti {WLUpL^)) 1 {tT{Up)tT{W)) 



pw 



MW)) 



MW)) 



(2.1) 



where Up = U^y{x) is the plaquette in the (/x, i^) plane. Note that the correlation 
function (|2.1j ) is sensitive to the field strength rather than to the square of the field 
strength [pT 



Pw 



> 



According to Eq.( |2.2| ) we define the color field strength tensor as: 



-F^u (x) 



pw[x) 



(2.2) 



(2.3) 



By varying the distance and the orientation of the plaquette Up with respect to the Wilson 
loop W, one can scan the color field distribution of the flux tube. 

We performed numerical simulations with Wilson action and periodic boundary con- 
ditions using an overrelaxed Metropolis algorithm. Our data refer to 16^, 20"^, and 24^ 
lattices. To evaluate the correlator Eq. ( p.l|) we used square Wilson loop Lw x L]y, with 
Lw = L/2 — 2 (L being the lattice size), and rectangular Wilson loops L/2 x L/4. 

In order to reduce the quantum fluctuations we adopted the controlled cooling algo- 
rithm 
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It is known |^0[ that by cooling in a smooth way equilibrium conflgurations, 
quantum fluctuations are reduced by a few order of magnitude, while the string tension 
survives and shows a plateau. We shall show below that the penetration length behaves 
in a similar way. 

For reader convenience let us, briefly, illustrate our cooling procedure. The lattice 
gauge conflgurations are cooled by replacing the matrix f/^(x) associated to each link 
/ = (x, fi) with a new matrix U'^{x) in such a way that the local contribution to the lattice 
action 



S{x) = 1 - -tT{U^{x)k{x)F{x)} 



(2.4) 



is minimized. F{x) = k{x)F{x) is the sum over the "U-staples" involving the link / and 
k{x) = ^ det ^F(x) j , so that F{x) G SU(2). In a "controlled" or "smooth" cooling step 
we have 



U^{x) 



V{x)U^ 



x] 



(2.5) 
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where V{x) is the SU(2) matrix which maximizes 

tT{V{x)Uf,{x)F{x)} (2.6) 
subjected to the following constraint on the SU(2) distance between U^{x) and U'^{x): 

\tT [{Ulix) - U'Jix)) {U,{x) - U'^ix))] < 6' . (2.7) 

We adopt 6 = 0.0354. A complete cooling sweep consists in the replacement Eq. (|275| ) 
at each lattice site. We do the above replacement vector-like according to the standard 
checkerboard order. 

The cooling technique allows us to disentangle the signal from the noise with a rel- 
atively small statistics. After discarding about 3000 sweeps to insure thermalization we 
collect measurements on configurations separated by 100 upgrades for 9 different values 
of (3 in the range 2.45 < [3 < 2.7 After cooling we obtained a good signal for on very 
small statistical samples (20 ^ 100 configurations). 

In Figure ^ we report our results for the field strength tensor F^y{xi,Xt), where the 
coordinates x/, Xt measure respectively the distance from the middle point between quark 
and antiquark (which corresponds to the center of the spatial side of the Wilson loop W 
in Eq. ( |2.1|) ) and the distance out of the plane defined by the Wilson loop. 

The entries in Fig. ^ refer to measurements of the field strength tensor taken in the 
middle of the flux tube {xi = 0) with 8 cooling steps ai (3 = 2.7 on the 24^ lattice, using 
a square Wilson loop W of size 10 x 10. Our results show that pw is sizeable when 
Up and W are in parallel planes. This corresponds to measure the component Ei of the 
chromoelectric field directed along the line joining the qq pair [E^ in Fig. |^). Moreover we 
see that Ei{xi, Xt) decreases rapidly in the transverse direction Xt- In Figure |^ we display 
the transverse distribution of the longitudinal chromoelectric field along the flux tube. 
The static color sources are at xi = +5 and xi = —4 (in lattice units). Figure ^ shows 
that the effects of the color sources on the chromoelectric fields extends over about three 
lattice spacings. Remarkably, far from the sources the longitudinal chromoelectric field 
is almost constant along the q — q line. Thus, the color field structure of the q — q tube 
which emerges from our results is quite simple: the flux tube is almost completely formed 
by the longitudinal chromoelectric field which is constant along the flux tube (if xi is not 
too close to the static color sources) and decreases rapidly in the transverse direction. 



B. Maximally Abelian projection 

In the 't Hooft formulation p the dual superconductor model is elaborated through the 
Abelian projection. The idea is that the Abelian projected gauge fields retain the long 
distance physics of the gauge system. In particular the physical quantities related to the 
confinement should be independent of the gauge fixing, and agree with those obtained 
in the full gauge system. This suggested us to investigate the Abelian projected 
correlator 

{trjW^U^)) l {tr{U^)tr{W^)) 
{tiiW^)) 2 (tr(Vr^)) ■ 
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The correlator is obtained from Eq.(|2.1|) with the substitution U^{x) — > U^{x). For 
instance the Abehan projected plaquette in the (/i, v) plane is 



= diag [exp {ze^^ix)) , exp (-zC(^))] • 



Obviously the Abelian projected quantities are commutating, so that we do not need 
the Schwinger lines in Eq.( |2.8| ). It is worthwhile to stress that is a gauge-dependent 
correlator. We performed measurement for 6 different values of (3 in the range 2.45 < 
(3 < 2.70 using the 16^ and 20^ lattices. In this case we find a good signal without 
cooling. Measurements are taken on a sample of 500 — 700 configurations each separated 
by 50 upgrades, after discarding 3000 sweeps to allow thermalization. The maximally 
Abelian gauge is fixed iteratively via the overrelaxation algorithm of Ref. [|ll|] with the 
overrelaxation parameter u = 1.7 (for further details see the Appendix). Remarkably 
enough, it turns out that the Abelian field strength tensor 

FU^) = ^pU^) (2.10) 

behaves like the gauge- invariant one defined by Eq. ( |2.3| ) . In Figure ^ we report our results 
for the field strength tensor F^y{xi,Xt) evaluated on maximally Abelian projected gauge 
configurations. The entries in Figure ^ refer to measurements done aX xi = -|-1 on a 
16^ lattice at /3 = 2.5 using a square Wilson loop of size 6 x 6 in Eq. ( p.8| ). Again 
we see that only the longitudinal chromoelectric field is sizeable. In Figure ^ we study 
the Xrdependence of the longitudinal Abelian chromoelectric field extracted using 6x6 



Wilson loop in Eq.(2J.) at = 2.5 on the 16^ lattice. Note that in the present case the 
static sources are at xi = +3 and xi = —2. The longitudinal Abelian chromoelectric 
field, likewise the non Abelian one, does not depend on the longitudinal coordinate xi, far 
from the static sources. It is worthwhile to observe that Fig. ^ suggests that the Abelian 
static sources are more localized than the non Abelian ones. This is in accordance with 
our previous observation that the maximally Abelian gauge fixing seems to reduce the 
fluctuations which are unimportant for the long distance physics. 

In the next Section we shall analyze our numerical data within the dual superconductor 
hypothesis. 



III. LONDON PENETRATION LENGTH 

A. SU(2) 

If the dual superconductor scenario holds, the transverse shape of the longitudinal chromo- 
electric field El should resemble the dual version of the Abrikosov vortex field distribution. 
Hence we expect that Ei{xt) can be fitted according to 

$ 

Ei{xt) = —p Ko{pxt) ,xt > (3.1) 

/TT 

where Kq is the modified Bessel function of order zero, $ is the external flux, and A = 
1/p is the London penetration length. Equation ( ^.1[ ) is valid if A 3> ^, ^ being the 
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coherence length (type-II superconductor). The length ^ measures the coherence of the 
magnetic monopole condensate (the dual version of the Cooper condensate). To determine 
the coherence length one should measure the correlation between the chromomagnetic 
monopoles. To do this one should construct a monopole creation operator. Unfortunately, 
thus far there is no a convincing proposal for the monopole operator. However, recently 
a promising proposal has been advanced in Ref. I^Tl- We shall return on this matter in 



Section |V] For the time being, because we are not able to determine the coherence length, 
we analyze our data far from the coherence region. To this end we try a fit with the 
transverse distribution ( |3.1| ) by discarding the points nearest to the flux tube {xt = 0). 

Let us discuss, firstly, the gauge invariant correlator Eg. ( |2.1|) . We fit Eg. to our 
data for Xt > 2 (in lattice units) obtaining x^/f ^ 1 ('^e used the MiNUlT code from 
the CERNLIB). In Figure ^ we show Ei{xt) measured in the middle of the flux tube 
together with the result of our fit. The fit results into the two parameters $ and fi. We 
have checked the stability of these parameter by fitting Eg. (|3.1|) to the data with the cuts 
Xt > a;™"^, X™" = 2, 3, 4, 5. In Table |l] we report the results of our stability analysis. We 
can see that within the statistical uncertainties the fit parameters are guite stable. So 
we are confident that our determination of the London penetration length is trustworthy. 
We ascertained, moreover, that the data obtained from the gauge-invariant correlator 
with cooled gauge configurations leads to a parameter /i which shows a plateau versus 
the number of cooling steps (see Fig. [^). This corroborates our expectation that the 
long range physics is unaffected by the cooling procedure. On the other hand. Figure || 
indicates that the overall normalization of the transverse distribution of the longitudinal 
chromoelectric field is affected by the cooling. In fact the parameter $ does not stay 
constant with the cooling. We feel that this is an indication that the flux $ is strongly 
affected by lattice artefact. This point will be thoroughly discussed below. 



In Figures ^ and we display the inverse of the penetration length (in units of 
Ajjg) and the external flux $ versus p. These data are obtained by fitting Eg.(|3.1|) to 
the data extracted from sguare Wilson loop (open points) and rectangular Wilson loops 
(full points). A few comments are in order. A look at Figure |^ shows that the inverse 
of the penetration length /x agrees within statistical fluctuations for both kinds of Wilson 
loops. However we see that for (3 > 2.65 the parameter /i arising from the rectangular 
Wilson loops seem to display sizeable finite volume effects. On the other hand we find 
that the parameter n extracted from the sguare Wilson loops displays finite volume effects 
for /3 > 2.7, in the case of the 24^ lattice. So in order to simulate in the range /3 > 2.7 we 
need lattices with L > 24. 

Figure | suggests that the ratio fi/Ajjg displays an approximate plateau in /3. Indeed 
we fitted the ratio with a constant and obtained 



8.96(31), x7/ = 2.11. (3.2) 



using sguare Wilson loops in Eg. (pTTI), and 



^ = 9.36(29), xV/ = 0.53. (3.3) 



^AIS 



for rectangular Wilson loops (discarding in the fit the points at j3 > 2.65). 

Eguations(^.2|) and(|3.3|) corroborates our previous observation on the consistency of 
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the penetration length. By fitting all the data we obtain 



-^ = 9.17(21), xV/ = l-48. (3.4) 

It is worthwhile to stress that our evidence for asymptotic scaling of the penetration 
length is only indicative. In general it should be much easier to check scaling rather than 
asymptotic scaling. We looked at the scaling of /i extracted from square Wilson loops 
with the square root of the string tension (we have used the string tension extracted from 
large Wilson loops). We found that there is approximate scaling of with ^/a for P > 2.5: 

-^ = 4.04(18), xV/=l-38. (3.5) 



So we see that our data on the penetration length are in agreement with the general 
expectation that scaling goes better than asymptotic scaling. On the other hand, the 
approximate evidence of asymptotic scaling is a natural consequence of the fact that the 



penetration length is a physical quantity related to the size D of the flux tube [1/ 



D ^ - . (3.6) 

As concern the parameter $, Figure |l^ shows that $ is rather insensitive to the shape 
of the Wilson loops used in Eq.( |2.lD (again the data from rectangular Wilson loops are 
affected by finite volume effects for (3 > 2.65). Moreover $ decreases rapidly by increasing 
P and seems to saturate to a value quite close to 1. We postpone the discussion of 
this behaviour until the comparison with the results obtained using Abelian projected 
configurations in the maximally Abelian gauge. 



B. Maximally Abelian projection 



Let us consider, now, the Abelian projected field strength tensor Eq. ( p.lO|) . As we saw, 
only the longitudinal Abelian chromoelectric field is sizeable. As in previous case we try 
to fit the data with the law: 



^ ,.2 



2tt 



(3.7) 



Again we find (see Fig. |TT]) that Eq. ( ^TTf ) reproduces quite well our data for Xj > 2 
ix^/f ~ !)• Table ^ we check the stability of the fit parameters. In Figure |T2| we 
display the ratio ^aI ^ws obtained by fitting Eq. (|3?7|) to the data in the case of square 
Wilson loops (open points) and rectangular Wilson loops (full points). Within the (rather 
large) statistical uncertainties the parameter agrees for the two different Wilson loops. 
Moreover the data suggest that the ratio iiaI ^tis '^^^^ depend on (3. Indeed we fit 
the ratio with a constant and find 



A 



MS 



using square Wilson loops in Eq. 

IJ'A 



A 



8.26(67) , 

I^^TSl) , and 
8.27(52) , 



XV/ = 0.41 



xV/ = l-87. 



(3.8) 



(3.9) 



MS 
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for rectangular Wilson loops. An overall fit of all the data gives 



-^ = 8.27(41), xV/ = l-05. (3.10) 

Note that Eq. ( p.2|) and Eqs. (|3.8| )- (|3rT0|) give consistent value for the ratio fi/Ajj-g. On the 
other hand the ratio fi/Aj^sj Eq. ( p.3|) , extracted from the gauge invariant correlator pw 
with rectangular Wilson loops is slightly higher than Eqs. ( |3.8| )-( |3rT0| ). Indeed, Eq. ( |3.3| ) 
and Eqs. ( p.8|) -( 3.10 ) are consistent within two standard deviations. We feel that this 



small discrepancy is due to the fact that the rectangular Wilson loops seem to be more 
sensitive to finite volume effects. For this reason we shall, henceforth, refer to the data 
extracted from square Wilson loops. 

In Fig. |13| we report the ratio fi/Aj^ and pa/ ^ws versus (3 obtained by the data 
corresponding to square Wilson loops. We can see that the London penetration length 
extracted from the gauge-invariant correlator Eq. ( |2.1| ) agrees with the one extracted from 
the Abelian projected correlator Eq. ( |2.8|) . In Figure |l3|we show also the result obtained 
by fitting together the data (for square Wilson loops): 

^ = 8.84(28), xV/ = l-44. (3.11) 



As a consequence we can safely affirm that the London penetration length is gauge in- 
variant. We feel that this results strongly supports the dual superconductor mechanism 
of confinement. 

As concern the parameter we find that, unlike the previous case, does not 
depend strongly on {3 (see Fig. |I4D. Moreover we see that $a is quite close to 1. It is 



worthwhile to discuss the physical interpretation of $. The total flux $t of the flux tube 
chromoelectric field is given by 

^T = j d'xtEiixt) , (3.12) 

where the integral extends over a plane transverse to the line joining the static color 
charges. As we have already discussed, the transverse distribution of the longitudinal 
chromoelectric field can be described by the law Eq. ( p.l| ) when Xt > 0. Obviously we 
cannot extend the validity of Eq. ( |3.1D up to — > 0. Indeed for — > we encounter a 
logarithmic divergence in Kq. On the other hand, Ei{xt) is finite in the coherence region 
Xt < ^. However, if A/^ > 1 we extimate that the extrapolation up to the origin introduces 
an overextimation of the integral (|3.12|) by less than 10%. So, inserting ( p.l|) into ( 3.12|) , 
we get: 

$T= /" d^XtEiixt)^^ . (3.13) 



Equations ( p.l2 ) and ( p.l3 ) tell us that the parameter $ measures the total flux if A/^ ^ 1. 



In U(l) it turns out that $ = 1, since that happens to be one unit of quantized electric 
flux [^. If the dynamics of the Abelian projected fields resembles the gauge fields of 
U(l), then we expect that $a — 1- Indeed we find (square Wilson loops) 

$A = 1.15(5), xV/ = 0.79. (3.14) 
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From the previous discussion it follows that Eq. ( p.l4|) seems to indicate that A/^ ~ 1. 
We would like to contrast Eq. ( 3.14] ) with the behaviour of In Figure we report 
and $ versus (3. The behavior of $ under the cooling (see Fig. ^ suggested that 
the external flux is strongly affected by lattice artefacts. Moreover, Figure |15] indicates 
that the lattice artefacts seem to disappear by increasing f3. Thus we are led to suspect 
that the external flux gets renormalized by irrelevant operators, whose effects are strongly 
suppressed in the maximally Abelian gauge. 



IV. STRING TENSION 

In the previous Section we have shown that the color fields of a static quark-antiquark 
pair are almost completely described by the longitudinal chromoelectric field. Moreover 
we showed that the longitudinal chromoelectric field is almost constant along the flux 
tube. This means that the long distance potential which feels the color charges is linear. 
Obviously the string tension is given by the energy stored into the flux tube per unit 
length. As a consequence we can write 



a 



^ I d'xtEf{xi,Xt) . (4.1) 



We stress that the string tension a defined by Eq. ( [4.1|) does not depend on xi as long as 
the longitudinal chromoelectric field is constant along the flux tube. As we have already 
discussed, working on a finite lattice results in the limitations xi = 0, ±1 (in lattice units) 
in the integrand in Eq. ( [4.1|) . Keeping these limitations in mind, from Equation we 
can obtain an explicit relation between the string tension and the parameters $ and /i. 
Indeed, if we extrapolate Eq. (|3.1| ) up to Xt = 0, by using 

X 

dxxK^{x) = -, (4.2) 



we get 



:/i . (4.3) 



The main uncertainty in Eq. (|4.3| ) comes out from the parameter $. As explained in 
Sects. II and III we computed the parameters $ and fi on SU(2) gauge configurations 
and on the maximally Abelian projected gauge configurations. In the latter case $^ ~ 1 
and independent of (3. On the other hand, for SU(2), $ > 1 and it approaches values 
very close to $^ by increasing (3. As we have already discussed, we feel that the external 
flux $ is strongly affected by lattice artefacts. We can try to get rid of these effects by 
assuming that in the limit /3 — oo 

$ ~ <1>^ ~ 1 . (4.4) 

In this way Eq. ( |4.3| ) becomes 

a ^ ^ . (4.5) 

Von 
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A striking consequence of Eq. ([4 .51) is that, due to /i ~ /i^, 

v^^v^A (4.6) 

within statistical uncertainties. 

In Figure |16| we report Eq. ([4.5|) in units of Ajjg versus aAj^. Fitting all together 



the data to a constant we get (square Wilson loops) 



1.76(6), xV/ = l-44. (4.7) 



The quoted error in Eq. (O) is purely statistic. However, one should keep in mind that 
our theoretical uncertainties in the extimation of the string tension ( |4.7| ) introduce a 
systematic error which can be of the order of ten per cent. Nevertheless, it is gratifying to 
see that our extimation of the string tension Eq. ( |4.7] ) is consistent with (star in Fig. |1^) 

^ 1.79(12) . (4.8) 



The value quoted in Eq. ( [4.8| ) has been obtained in Ref. [^] by the linear asymptotic 
extrapolation of the string tension data extracted from Wilson loops on lattices larger 
than ours. 



V. CONCLUSIONS 

Let us conclude by stressing the main results of this paper. We investigated the color 
field strength tensor of the q — q flux tube by means of the connected correlators (|2.1| ) 
(full SU(2)) and ( |27^ ) (maximally Abelian gauge). 

The main advantage of using the connected correlator ( p.l|) and ( p.8|) resides in the 
fact that the connected correlators are sensitive to the field strength rather than to the 
square of the field strength. As a consequence we are able to detect a sizeable signal even 
with relatively low statistics. It turns out that the flux tube color fields is composed by 
the chromoelectric component parallel to the line joining the static charges. Moreover 
the longitudinal chromoelectric field is almost constant far from the color sources, and 
it decreases rapidly in the directions transverse to the line connecting the charges. As a 
matter of fact we found that the transverse distribution of the longitudinal chromoelectric 
field behaves in accord with the dual Meissner effect. This allows us to determine the 
London penetration length. We checked that the penetration length is a physical gauge 
invariant quantity. A remarkable consequence of our findings is that the long range 
properties of the SU(2) confining vacuum can be described by an effective Abelian theory. 



In addition, after fixing the gauge with the constraints (|1.8|) - (|1.10|) , it seems that the 
degrees of freedom which are not relevant to the confinement get suppressed. 

Finally, we put out a very simple relation between the string tension and the pene- 
tration length which gives an extimate of ^/a quite close to the extrapolated continuum 
limit available in the literature. 

In conclusion we would like to stress that the most urgent problem to be addressed 
in the future studies is the reliable extimation of the coherence length C,- The results in 
Sect. Ill give an indirect and, admittedly, very weak indication that A/.^ ~ 1. As we 
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have already discussed, the coherence length is determined by the monopole condensate, 
the order parameter for the confinement. Recently two different groups |^,^ give 
an extimation of the coherence length. These authors calculate the electric fiux and 
magnetic monopole current distribution in the presence of a static quark- ant iquark pair 
for SU(2) lattice gauge theory in the maximally Abelian gauge. The magnetic monopoles 
are identified using the DeGrand-Toussaint construction. By using a dual form of the 
Ginzburg-Landau theory |^|, which allows the magnitude of the monopole condensate 



density to vary in space, they fitted the data and obtain A/a and ^/a. They found that 
the coherence length is comparable to the penetration length. Even though we feel that 
the approach of Refs. [p4| , p5[ is interesting, we would like to observe that it relies heavily 
on the definition of the magnetic monopole current. As a matter of fact, in Ref. ||28| it was 
pointed out that the DeGrand-Toussaint definition of the monopole density is plagued by 
lattice artefacts, which are, however, less severe in maximally Abelian gauge. So that the 
DeGrand-Toussaint monopole density is not an order parameter for confinement. Thus 
the approach of Refs. [|^, is plagued by the ambiguities related to the definition of the 
monopole current. On the other hand, in our approach we work outside the coherence 
region, so that we feel that our results do not manifest the above mentioned problem. To 
clarify this point, it should be of great help the study of the distribution of color fields 
in the presence of a static quark-antiquark pair in the framework of the dual Ginzburg- 
Landau model with the magnetic monopole current constructed by means of the monopole 
creation operator proposed in Ref. . 



APPENDIX 

In this Appendix we give more details on the algorithm used to fix the maximally Abelian 
gauge. On the lattice, the maximally Abelian gauge is obtained by maximizing the lattice 
functional 

Ri = J2l^^hU,ix)a,UjXx)] (Al) 

X,fi 

over all SU(2) gauge transformations 

U,{x) U,{x) = gix)U,{x)g^{x) , (A2) 

where (7 (x) G SU(2). Under an arbitrary gauge transformation the variation of the lattice 
functional Eq. ( [Al| ) is 

ARi{x) = hi [g\x)(T3g{x)X{x)] - hi [a^Xix)] , (A3) 

where 

X(x) = J2 [U,ix)a,Ulix) + Ulix - fi)a,U,{x - /i)] (A4) 

belongs to the SU(2) algebra. If we have locally maximized the lattice functional ( |AT| ) 
with respect to an arbitrary gauge transformation, then we have 

ARi{x) = . (A5) 
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From Eq. ( [A3| ) it follows that 

X{x) = g{x)X{x)g\x) , (A6) 

i.e. X{x) must be diagonal. So that maximizing Ri{x) is equivalent to diagonalizing the 
hermitian matrix X{x). Note that maximization of Ri{x) at the given lattice site x is 
accomplished by a gauge transformation g{x), which, in turn, affects the value of the 
local operator X{x) at the nearest neighbours. Therefore the maximization of the lattice 
functional (|A1|) can be achieved only by an iterative procedure over the whole lattice. 
In equivalent manner one can find g{x) as the matrix which diagonalizes X{x) or as the 
matrix which maximizes Ri{x). 

To obtain explicitly the gauge element g{x) which maximizes Ri{x), let us write Ri{x) 

as 

Rii^) = ^tr [(T3g\x)(T3g{x)X{x)a3] 

2 ^ (A7) 
= k{x)-ti [a3g\x)a3g{x)V{x)] 



where 



Vix) = , k{x) = Vdet {X{x)as) (A8) 



ensuring that V{x) is an element of SU(2). As one can easily recognize from Eq. 

V{x) = vo{x) + i{vi{x)cri + V2{x)a2) ■ (A9) 

Note that the term proportional to (T3 is absent. Now, we observe that, if we consider 
g{x) = u{x)g{x) (with u{x) = Uo{x) +iu3{x)(r3) instead of g{x), then Eq. ( |A7| ) is invariant. 
So we can assume without loss of generality that in Eq. ( |A7|) g{x) = go{x) + i{gi{x)ai + 
g2{x)(J2). As a consequence, Eq. (A7) is maximized when 



Vo{x +1 
go{x) = ±y , 

gi{x) = ~^ (AlO) 
92[x) - 



2go{x) 



Since maximization of ( |Al| ) results in an iterative procedure, we must have at disposal a 
convergence criterion. To have a measure of the goodness of gauge fixing we consider the 
average size of the non-diagonal matrix elements of X over the whole lattice: 

^'^'r)=7iEo^ii'+i^2r] (All) 



L4 

X 

where X = Xiai + Xia2 + Xia^. We stop the iterations when 

\X'"^\^)<D (A12) 
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where D is some (small) positive number. In our simulation we used D = 10~^. 

In order to accelerate the convergence of the algorithm we adopted the overrelax- 

Once we have found the matrix g{x) which 



ation method 
maximizes (\K 



suggested in Ref. |]TT 



we make the following substitution 



9o 



(A13) 



where the overrelaxation parameter u varies in the interval 1 < uj < 2. The exponentiation 
in Eq. ([A13|) is obtained through the following representation for an element u G SU(2): 



u = cos(^) + 2 (a ■ f) sin(^) , 



(A14) 



where f = r/|?^. 

In our Monte Carlo runs we used u = 1.7. However, we would like to stress that it 
exists [|11| an optimal overrelaxation parameter Uc- Moreover for large lattice size L it is 
believed that 



(A15) 



where the constant c is problem dependent. As a matter of fact it turns out that a better 



convergence can be obtained using values of uj close to 1.9 (see Figure |T^. Indeed we 
obtained Uc — 1.92 for L = 16. Inserting this value into Eq. (|A15|) we find c ~ 0.7. It is 
remarkable that our value for the constant c agrees with the one relevant to the Landau 



gauge fixing [jTT 
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(3 = 2.7, using Wilson loops of size 10 x 10 in Eq. ( pi]) . 
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Figure 6: The London fit ( p.l|) to the data for the longitudinal chromoelectric field. 
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Figure 7: The inverse of the penetration length /x versus the number of the coohng steps 
obtained by fitting the transverse profile of the longitudinal chromoelectric field at xi — 
(8x8 Wilson loop). 
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Figure 8: The parameter $ in Eq. obtained by fitting the transverse profile of the 
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points refer to 8 x 8 Wilson loop, full points to 10 x 5 Wilson loop. 
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Figure 11: London fit ( p.7| ) to the data for the Abehan longitudinal chromoelectric field 
at X; = for square Wilson loop. 
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Figure 12: /x^/ A versus /3. Full points correspond to rectangular Wilson loops, open 
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Figure 16: String tension (in units of A;^) evaluated through Eq. ([4.1[). Star refers to the 
value given in Ref. |2^. Symbols as in Fig. |15|. For figure readability not all the available 
data are displayed. 
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GAUGE FIXING SWEEPS 

Figure 17: Efficacy of gauge fixing defined by Eq. ( [All| ) as a function of the overrelaxation 
parameter uj for the L = 16 lattice. The case uj = uj* corresponds to alternate uj = 1.0 
with UJ = 2.0 in the gauge fixing sweeps. 
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fit parameters stability 
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5 
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Table 1: Fit parameters in Eq. ( p.l|) versus of the number of discarded points in the Xt 
direction. 
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Table 2: Fit parameters in Eq. ( p.7| ) versus of the number of discarded points in the Xt 
direction. 
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